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ABSTRACT 


The formalism of Supersymmetric Quantum Mechanics provides us the eigenfunctions to 
be used in the variational method to obtain the eigenvalues for the Hulthen Potential. 
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1 Introduction 

The Hulthen Potential, in atomic units, is given by: 


V H (r) 


Se~ Sr 
1 - e~ Sr ’ 


( 1 ) 


where 5 is the screening parameter. This potential has been used in several branches of 
Physics, (see [jlj and references therein). 

The Schroedinger equation for this potential can be solved in closed form for s waves, 
0 . However, for l ^ 0 it is necessary to use numerical methods, [|l|] . 

Supersymmetric Quantum Mechanics (SQM) has been used to solve Schroedinger equa¬ 
tion of solvable potentials, [[|, partially solvable ones, [|J, in the WKB-approximation, j[|. 
and it has also been applied in variational method, as recently suggested, @ , [0] • 

The supersymmetric formalism has already been used to study some aspects of the 
Hulthen potential, ||, [hj. Here, the exact analytical solution for this potential is reob¬ 
tained, for l = 0 to show the consistency of the method. When I ^ 0 we interpret that the 
supersymmetry is “broken” by the potential barrier terms. Nonetheless, the supersymmetry 
gives us eigenfunctions that allow us to compute the eigenvalues of the variational method. 
The eigenfunctions for 2 p and 3d states are evaluated for some values of the parameter delta. 
Our results are compared with direct numerical integration data, JI|. 


2 The Hulthen Potential with l = 0 

The Hamiltonian for the Hulthen potential (l = 0) can be written as: 

H - - 6e ~ Sr 

1 2 dr 2 1 — e~ 5r 

From the Hamiltonian hierarchy we can obtain its eigenfunctions and eigenvalues, 
In this way, first we factorise the Hamiltonian (||): 

Hi - E ( o ] = a+a^ 

with creation and anihilation operators defined by 


( 2 ) 

IT 


( 3 ) 


( 4 ) 


and we determine the ground state eigenvalue E^\ the superpotential W\{r) and the ground 
state eigenfunction 

r fo Wl ( r ') dr '. (5) 


To 1} oce h 


From the operators af and aq , we write H 2 ,the supersymmetric partner of Hi as 


H 2 - E^ = eqa+ 


( 6 ) 
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and factorise H 2 in an analogous way to (Q), 


a 2 a 2 = H 2 — E, 


( 2 ) 


( 7 ) 


with a 2 defined similarly to (ffj). We then determine E 0 , W 2 (r) and Tq . Interchanging 
the operators a 2 and a 2 in (0) we obtain the supersymmetric partner of H 2 . Repeating 
this process n times we obtain the eigenvalues and eigenfunctions for the n-th ground state 
Hamiltonian. The superalgebra allows us to relate these results with the original Hamiltonian 
by the relations: 

E m = £<»+■>, *M( r ) = 0 +a+..oi*[," +1) (r). 


In our case, using the Hulthen potential, the n-th superpotential is given by 

nSe~ 5r 


WJr = 


o—5r 


1 n 
u~2 6 


that corresponds to the n-th member of the Hamiltonian hierarchy: 


Vn(r) = Wl(r) 


d ... . . n(n - 

= “ 2(1 


1)5 2 


0 — 28r 


[n( 1 — n)S + 2 ]8e 


—Sr 


The energy-eigenvalue and eigenfunction, eq.( 

= If. 


- e~ Sr ) 2 
|), are 
n r 1 

2 > s + ny 


2(1 


o—5r 


Tj, n) (r) = (1 - e~ 5r ) n e- [ ^ 5]r . 

We can verify that the energy-eigenvalues (|lTD are the same given in reference 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 

( 12 ) 


3 The Hulthen Potential with l / 0 

The Hamiltonian for the Hulthen potential when l ^ 0 is written as: 

1 cP Se~ 3r l (l + 1) 


H = - 


2 dr 2 


0 —Sr 


2 r 2 


(13) 


The potential barrier term prevent us to build the superfamily as in the l = 0 case, 
since the potential is not exactly solvable. However, several numerical approaches have been 
used in order to evaluate the spectra of energy-eigenvalues and eigenfunctions. In particular, 
Greene and Aldrich, [12[], suggested an Hulthen effective potential in which the eigenfunctions 
are used as variational trial wave functions. 

Based in our results for the case l = 0, we introduce a new effective potential whose 
functional form is suggested by eq. ( 0 ), 


V ell M = 


5e 


—5r 


2 2<5r 


1 — e 


—5r 


1(1 + 1) S 2 e 
2 (1 - e-^) 2 ' 


(14) 
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We note that, for small values of 8, the second term of (13) gives us a potential barrier 
term of (|T3|) in first approximation. This potential differs from the one used in 0 by 
the exponential numerator which is a quadratic term. The advantage of using the effective 
potential (|T4D is the fact that we can vary its parameters without changing its functional 
form. This is reinforced by the constructive method of determining wave functions based on 
sup ersy mmet ry. 

As the effective potential given by (|b|) has the same functional form as (0). we can 
solve the Schroedinger problem by the factorisation method of SQM and find the whole 
super family. The superpotential of the n-th member of the family is given by 

s,—5r 


W n (r) = B n 


1 


o—5r 


+ C n 


where 


B n — — + 


<5 2 +4 Bn.^Bn^-5), C n = 


B n -\{5 — 2C n -\) — B n 8 


Bi — --(1 + \Jl + 4/(7 + 1)), 


and the energy and wave functions are given by 

5 

2 

The energy eigenvalues and wave functions given by 

1 
2 


(||), where 
n _ 8 B\ + 2 


2 Bn 


B n = 0 . 


|) are 


rp(n+l) _ x ^2 

-C' n — — 


n +1 


= (1 - e-") 


Cjr 


(15) 

(16) 

(17) 

(18) 
(19) 


Notice that B's and C' s depend on l. 

Thus, the spectrum of the effective potential (|T^) has quantum number n — 0,1,... that 
shows us which member of the superfamily is going to be used. For instance, n — 0 and 
l = 1 corresponds to state 2p\ n = 1 and l = 1 corresponds to state 3d and so on. Fixing 
l — 0 we recover the original Hulthen potential, as expected. 

We now look at the Hulthen potential (|13|) and solve the Schroedinger problem by the 
variational method. We start from V e ff given by ( ]T4|) whose eigenfunctions are given by (§) 
and (|l9l). For the state 2 p we use the first member of the superfamily (Z = 1) as our trial 
wave function, changing 8 by the variational parameter p, i.e., 


= si, 1 ’ 






( 20 ) 


The energy is obtained by minimisation with respect to p. Thus, the equation to be 
minimised is 

5e~ Sr , l{l+ 1) 1 


E,= 


r°° vD _ 

JO /^ \ / [ 2 dr ^ 1—e ^ 


+ 


2 r 2 


-]T M (r)dr 


/ 0 ~*„(r)Vr ' (21) 

The second derivative in (^1 ]) can be evaluated analitically. However, the integration has 
to be carried out numerically. Our explicit values for the 2 p, (l = 1) and 3d, (/ = 2) energy 
states for some values of the parameter 8 are listed bellow in Table 1. They are shown 
together with direct numerical integration data. 
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Table 1. Energy eigenvalues as a function of the screening parameter for the states 2 p 
and 3d, [eq. (|2T|) 1. Comparison is made with numerical data of Ref. [[L]]. 


State 

Delta 

Variational result 

Numerical Integration 

2p 

0.025 

-0.112760 

-0.1127605 


0.050 

-0.101042 

-0.1010425 


0.075 

-0.089845 

-0.898478 


0.100 

-0.079170 

-0.0791794 


0.150 

-0.059495 

-0.0594415 


0.200 

-0.041792 

-0.0418860 

3p 

0.025 

-0.043601 

-0.0437069 


0.050 

-0.032748 

-0.0331645 


0.075 

-0.023010 

-0.0239397 


0.100 

-0.014433 

-0.0160537 
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4 Conclusions 


We have obtained from the formalism of SQM the exact analytical eigenfunction and energy 
eigenvalue for the Hulthen potential for l = 0. When l ^ 0, we used an effective potential 
suggested by the case l = 0, to give us a variational trial wave function. The energies for the 
2 p and 3 d states were obtained for some values of the parameter 8. 

We have used the same variational trial wave function for the 2 p and 3d states, eq.(JT^). 
The only change was the angular moment number l. However, as the suggested potential 
has the same functional form, we could determine all the eigenfunctions and therefore we 
could use the most appropriate in each case, (see 


13 


We note that better results have been obtained for the 2 p state, for small values of 5. 
This is expected since for small values of 8 the effective potential (14) becomes closer to 
the original Hulthen potential (|I3|), and for l small the contribution of this angular moment 
term in the potential is also small. The advantage of this process is that only the form of 
the potential (0) is important to obtain the eigenfunctions; the multiple parameters can be 
changed conveniently in order to obtain better results. 

The algorithm used here is very powerful because it gives us a systematic method to look 
for variational trial wave functions based on the functional form of the effective potential in 
the factorisation method of SQM. 
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